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Abstract 

The distribution function for the intensity of radiation propagating in a ran- 
dom medium is analyzed for arbitrary multiphcity of scattering (for arbitrary 
relation between the distance of propagation and mean free path), including 
as limiting cases purely ballistic and purely diffusive regimes. The calcula- 
tion is fulfilled in the framework of diagram techniques and is based on the 
topological analysis of diagram. 

When electromagnetic wave propagates in a random medium and undergoes multiple 
scattering the density of energy in a given point (intensity) I =< > is a, strongly 
fluctuating quantity. In a diffusive regime these fluctuations are traditionally described by 
Rayleigh statisticsB. The distribution function Pr{I) is a simple negative exponential: 

PH(/) = ^exp(-^). (1) 

which corresponds to the equation for the moments: 

< /" >= n! < J . (2) 

In Eqs. (P and (^) the averaging is with respect to all macroscopically equivalent but 
microscopically different scatterers configurations. By diffusive we mean the regime when it 
is possible to discard the "coherent" component of the field, that is 



< ^ >= 0. (3) 

This assumption is true for the case r/£ ^ 1, where r is a distance from the source to the 
observation point, and £ is a mean free path (for the sake of definiteness we consider in 
this paper only the case of point source in an infinite medium). The aim of the present 
pubhcation is to get the distribution function for arbitrary relation between r and when 
generally speaking the "coherent" component of the field should also be taken into account. 

Here some clarification should be made. In our previous publicationJl we presented 
systematic derivation of the intensity distribution function in a diffusive regime and have 
shown that the Rayleigh statistics is correct only for // < I ><^ kl, where k is the radiation 
wave vector. For larger values of intensity the distribution function differs drastically from a 
simple exponential and the asymptotic behavior is a stretched exponential. In this paper we 
consider the case of M is so large that this effect can be ignored (or alternatively we are not 
interested in the distribution function tail). But if in the above mentioned publications we 
restricted ourselves by the case r/£ ^ 1, now we are interested in the case of the arbitrary 
value of the parameter r/C.. 

To achieve this aim we apply the same kind of reasoning based on the topological analysis 
of the diagrams we used previouslyi. To get the idea of it let us remember how Rayleigh 
statistics is obtained in the framework of the traditional diagram technique!!. The point 
source is situated at the origin; the intensity is measured in the point r. In the diagram- 
matic representation < / > is given by the diagrams with a pair of wave propagators 
and G^^, summed with respect to all possible interactions with the scatterers. The n-th 
moment < /" > is given by the set of diagrams with n propagators and n propagators 
Gq^. For us is important the following property of diagrams: if the diagram consists from 
several disconnected parts then the contribution of that diagram is equal to the product 
of contributions of disconnected parts. It means, that if we consider for < /" > only the 
diagrams consisting of n disconnected parts, each part being the set of diagrams with a pair 
of propagators (advanced and retarded), we immediately obtain Eq. (g). The multipher n! 
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which appears in the n-th moment is of combinatorial origin; it is simply the number of 
possible pairings between propagators. 

So to get Rayleigh statistics we have ignored the connected diagrams with two or more 
pairs of propagators. In Ref. we have formulated a perturbation theory which systemat- 
ically takes such diagrams into account, ant it is such diagrams which lead to deviations 
from negative exponential statistics for large values of intensity. But we also ignored the 
diagrams which contain isolated (dressed) propagators. The ground for this is that each 
isolated propagator gives contribution of the order of exp(— r/£), and in the diffusive regime 
one can ignore such diagrams. Here we should take them into account. That is for < /" > 
we consider only diagrams which consists of isolated propagators and pairs of propagators 
(advanced and retarded one). The sum of all diagrams can be written down in the following 
way: 

n 

< /" >= ^ P(n, m)A"^5"-™, (4) 

m=0 

where A = \< E >f, B =< I > —\ < E > \'^, and P{n,m) is the coefficient of purely 
combinatorial origin: 

P(n,m)= -. (5) 

The distribution function P{I) is connected with its moments by usual way: 

p(/) = r exp(.e/) E <r> ^. (6) 

Substituting the Eq. (|) into the Eq. (H) and changing the order of summation,which is 
possible due to convergence of the series, we get: 

P{I) = / exp(^e/) E ^ ,d E {-liBY-"^ ^. (7) 

Taking into account that for n > m 

x"-'"= (8) 



(n-m)! dX' 



we obtain: 



(m!)2 „^(n-m)!^ ^ 1 + z^i? 1^ 1 + j 



So finally Eq. can be written in the form: 



^ ^ 7-00 1 - a) 2tx ^ ' 

(the intensity is measured in the units of < / > and parameter a = A/ {A + B) = \ < E > 
P/ < / >, is introduced). It is obvious that for a = (purely diffusive transport), the 
Eq. (|10D coincides with the Eq. (|l]), so this parameter gives the degree of deviation from 
Rayleigh statistics. In the opposite limiting case a = 1 the Eq. (p!0D gives P (/) = 5 (/ — 1). 
This limiting case means that there are no fluctuations, which can occur only when there is 
no scattering at all, in other words in the ballistic regime. 

For intermediate value of a the distribution function can easily be calculated numerically. 
This behavior is represented on Fig. 1. We see that even for a <^ 1 (almost diffusive 
transport) taking into account of the "coherent" component of the fleld drastically changes 
the distribution function for small /; in particular for J = we get: P = 1/2. Also interesting 
are the predicted oscillations of the distribution function. 
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FIGURES 

FIG. 1. Distribution function for: (1) a = .1 (solid line), (2) a = .5 (dashed line), (3) a = .9 
(dot-dashed line). 
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